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Abstract 

We investigate ideal quantum gases in D-dimensional space and confined in a 
generic external potential by using the semiclassical approximation. In par- 
ticular, we derive density of states, density profiles and critical temperatures 
for Fermions and Bosons trapped in isotropic power-law potentials. Form 
such results, one can easily obtain those of quantum gases in a rigid box and 
in a harmonic trap. Finally, we show that the Bose-Einstein condensation 
can set up in a confining power-law potential if and only ii D /2 + D /n > 1, 
where D is the space dimension and n is the power-law exponent. 

PACS numbers: 05.30.-d, 05.30.Fk, 05.30.Jp 



1 Introduction 



For dilute alkali-metal atoms in magnetic or magneto-optical traps at very 
low temperatures, the Bose- Einstein condensation has been achieved^ in 1995 
and the Fermi quantum degeneracy^ in 1999. These results have renewed the 
theoretical investigation on Bose and Fermi gases. 

In the experiments with Bosons, the system is weakly-interacting and the 
thermodynamical properties depend on the s-wave scattering length (for a 
review see Ref. 3). Nevertheless, by using Feshback resonances, it is now 
possible to modify and also switch-off the atom-atom interaction.^ In the 
case of Fermions, the s-wave scattering between atoms in the same hyperfine 
state is inhibited due the Pauli principle. It follows that at low temperature 
the dilute Fermi gas, in a fixed hyperfine state, is practically ideal. ^ 

In previous papers we analyzed ground-state and vortex properties of 
Bose condensates in different external potentials: harmonic potential, 
toroidal potential^^ and double- well potential. Recently, we have also af- 
forded the study of the thermodynamics of interacting Bose gases in harmonic 
potential. ^^'^^ 

In this paper, we investigate the thermal properties of both Bose and 
Fermi ideal gases in a generic confining external potential. All the calcula- 
tions are performed by assuming a D-dimensional space. Such an assumption 
is motivated by esthetic criteria but also by recent experiments with degener- 
ate gases in systems with reduced or fractal dimension.^ We analyze in detail 
the isotropic power-law potential, from which one easily deduces the results 
of a rigid box and a harmonic trap. 
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2 Confined Ideal Fermi and Bose Gases 



Let us consider a confined quantum gas of non-interacting identical Fermions 
(or Bosons) in D-dimensional space. In the grand canonical ensemble of 
equilibrium statistical mechanics/^ the average number No. of particles in 
the single-particle state \a) with energy is given by 

^ g/3(ea-/.) J- X ' 

where the sign + (— ) is for Fermions (Bosons), is the chemical potential and 
/3 = l/{kT) with k the Boltzmann constant and T the absolute temperature. 
In general, given the single-particle function N^, the average total number 
N of particles of the system reads 

N = J2Na. (2) 

a 

This condition fixes the chemical potential /i. Thus, // is a function of /3 and 
N. In the case of Fermions, fj, has no limitations and at zero temperature /j, 
is called Fermi energy Ep. From the Fermi energy Ep one immediately ob- 
tains the Fermi temperature Tp — Ep/k. Below the Fermi temperature, the 
Fermions begin to fill the lowest available single-particle states in accordance 
with the Pauli exclusion Principle: one has the Fermi quantum degeneracy. 
In the case of Bosons, n cannot be higher than the lowest single-particle en- 
ergy level eo, i.e., it must be < eo- When — > eo the function Nq diverges 
and consequently also N diverges. The physical meaning is that the low- 
est single-particle state becomes macroscopically occupied and one has the 
so-called Bose- Einstein condensation (BEC).^^ It is a standard procedure to 
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calculate the condensed fraction Nq/N and also the BEC transition temper- 
ature Tg by studying the non divergent quantity N — Nq at fj, — as a 
function of the temperature.^'^^ 

In the semiclassical limit, the D-dimensional system is described by a 
continuum of statcs-'^^^^'' and. instead of e^, one uses the classical single- 
particle phase-space energy e(r, p), where r = (ri, ru) is the position 
vector and p = (pi, ■■■,Pd) is the linear momentum vector. In this way one 
obtains from Eq. (1) the single-particle phase-space distribution 

Note that the accuracy of the semiclassical approximation is expected to be 
good if the number of particles is large and the energy level spacing is smaller 
then kT}^'^"^ Because of the Heisenberg principle, the quantum elementary 
volume of the single-particle 2D-dimensional phase-space is given by {27ih)^, 
where h is the Planck constant. It follows that the average number N of 
particles in the D-dimensional space can be written as 



N 

where 



d^p 

is the spatial distribution, and 

is the momentum distribution. It is important to observe that the total 
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number N of particle can also be written as 




1 



(7) 



g/3(e-/.) ± I 



where p(e) is the density of states. It can be obtained from the semiclassical 
formula 



where 5{x) is the Dirac delta function. 

In the case of Fermions, at zero temperature, i.e., in the limit /3 — > oo 
where /j, Ep (the Fermi energy), the phase-space distribution (3) becomes 



where 6(a;) is the Heaviside step function. ^^'^^ 

In the case of Bosons, below the BEC transition temperature Tb, the 
equation (3) describes only the non-condensed thermal cloud. Thus, the semi- 
classical quantization renormalizes the exact Bose distribution (1) that is di- 
vergent. Because there is not an unique way to translate wave-functions into 
a phase-space distribution,^^ one cannot introduce an exact single-particle 
phase-space distribution for the Bose condensate. Nevertheless, the exact 
spatial distribution of the Bose condensate is no(r) = |\E'(r)p, where \E'(r) is 
called order parameter or macroscopic wave-function of the condensate, nor- 
malized to the number Nq of condensed Bosons. For an ideal Bose gas, the 
function ^(r) is simply the eigenfunction of the lowest single-particle state 
of the system. In this paper we do not study the density profiles of the Bose 
condensed fraction because, for a non-interacting gas, their shape is not tem- 
perature dependent: only their normalization is a function of temperature. 





n(r,p) = e(£;F-e(r,p)) , 



(9) 
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Actually, to calculate the BEC transition temperature Tb and the condensed 
fraction Nq/N it is sufficient to study the non-condensed fraction (thermal 
cloud). (For a recent discussion of the properties of the Bose condensate, see 
Ref. 3 and also Ref. 5-14). 



3 Gases in External Potential 

Let us consider the ideal Fermi (Bose) gas in a confining external potential 
C/(r) and in a D-dimensional space. The classical single-particle energy is 
defined as 

<r,p) = |^ + t/(r), (10) 

where p^/ (2m) is the kinetic energy and m is the mass of the particle. The 
D-dimensional vectors r = [xi,...,X]j) and p = (pi,...,Pd) are respectively 
the position and momentum of the particle. 

First, we note that, by using Eq. (8) and Eq. (10), the semiclassical 
density of states can be written as 

2 

where T{n) is the factorial function. Then, we introduce the Fermi and Bose 
functions. ^^'^^ 



DEFINITION 1. The Fermi function is given by 
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and the Bose function is 

1 , ze-yy"-'^ 

9n{z) = / dy- — , 

L[n) Jo 1 — ze y 

where r(n) is the factorial function. 

REMARK 1. The Fermi and Bose functions are connected by the relation 
fn{z) = —gn{—z). For \z\ < 1 the two functions can be written as 

oo 

1=1 * 

oo i 

i=l * 

One finds also that QniX) = C{''^)> where C{n) is the Riemann (-function}^'^^ 

Now wc state two theorems about ideal Fermi and Bose gases in external 
potential. Remind that we work in the semi classical limit. 

THEOREM 1.1 For an ideal Fermi gas in an external potential U{r) and 
D-dimensional space, the finite temperature spatial distribution is given by 

where A = (27r/i^/3/m)^/^ is the thermal length and ji is the chemical potential. 
The zero temperature spatial distribution is 

where Ep is the Fermi energy. 
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Proof. One finds the spatial distribution from tlie Eq. (5) by integrating 
over momenta the Eq. (3) (with the sign +). In particular, one has 

d^p 1 



n(r) 



2 



DtI^ ^i?-lg/3(/.-C/(r))g-/3(iL-) 

dp 



where Dn^ /r(-^ + 1) is the volume of the D-dimensional unit sphere. Then, 
with the position y"^ = and using the Fermi function fD{z) with z = 
Ql3(p.-u{r)) ^ one gets the finite temperature spatial distribution. Finally, one 
obtains the zero-temperature result by observing that 



n r 



2T7T^n {^F - U(r)) dpp^ 



-^c/ yrjF - 'J yi-)) 1^ u,p ^ 

1 DttT „ _ 1 / . X D 



-e{EF-U{r))-iJ2m{Ep-U{r) 



(27r^2)^r(f + 1) ' ' ^ " D 
where the spatial distribution is taken from the Eq. (8). □ 

In the same way, but using the sign — in Eq. (3) and the Bose function 
go [z) with z — e^'^^"^^^^^ , one can easily prove also the following theorem. 

THEOREM 1.2 For an ideal Bose gas in an external potential U{r), the 
finite temperature non-condensed spatial distribution is given by 
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where A = (27r^^/3/m)^/^ is the thermal length and ji is the chemical potential. 
□ 

These two theorems are the generahzation, of the formulas for ideal ho- 
mogenous Fermi and Bose gases in a box of volume V (for D — ?, see Ref. 
15). They show that, in the semiclassical limit, the non-homogenous results 
are obtained with the substitution ji fi — U{r), also called local density ap- 
proximation. In particular, with U{r) = 0, from the previous theorems, one 
obtains the Fermi temperature Tp and the Bose temperature Tb for quan- 
tum gases in a rigid box, by imposing the normalization condition (4). The 
results are 



where n = N/V is the homogenous density of particles (again, for D = 3 see 



In general, to find the momentum distribution, the Fermi temperature 
and the Bose temperature, it is necessary to specify the external potential. 
In many experiments with alkali-metal atoms, the external trap can be accu- 
rately modelled by a harmonic potential.^ More generally, one can consider 
power-law potentials, which are important for studying the effects of adia- 
batic changes in the trap. The density of states of a quantum gas in the 
power-law potential U{v) — A can be calculated from Eq. (11) and reads 




and 




Ref. 15). 




(12) 
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We can now state two theorems about ideal Fermi and Bose gases in 
isotropic power-law potentials. 



THEOREM 2.1 Let us consider an ideal Fermi gas in a power-law isotropic 
potential U{r) — A r" with r = |r| = {^iLi^^Y^^ . The finite temperature 
momentum distribution is given by 

"(P'-(2ftv^)"r(f + i) [lu) )■ 

The zero temperature momentum distribution is 



{2hy/7r)^ \AJ \ 2m J \ 2m) ' 
The Fermi energy Ep and the Fermi temperature Tp are given by 



Ep = kTp 



'H?!)*.4vIi|±Ilrf^ + ^ + i|iv 

m + \2 n 



1 

D.D 



where N is the number of Fermions in the gas. 

Proof. One finds the finite temperature momentum distribution from the 
Eq. (6) and by integrating over space coordinates the Eq. (3) (with the sign 
+). In particular, one has 

d^r 1 



n 



dr- 



where again Dtt^ /V{^ + 1) is the volume of the D-dimensional unit sphere. 
Setting = (3Ar'^ and using the definition of Fermi function /d (2;) with z — 
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e^^^^~2^\ one finds the first formula of the theorem. The zero-temperature 
results are obtained by observing that 



1 -ef^.-^llr = 



(2rf)T(f + l) V ' 2m, 



(27r;i^)^ r(f + 1) V 2mj D\ ^ " 2m' 
where the momentum distribution is taken from the Eq. (9). The Fermi 

energy Ep and the Fermi temperature Tp are found from the normalization 

condition of the zero-temperature momentum distribution. Namely, one finds 

J ^{2h^yn)^\Aj \ 2m J \ 2m J 

2 

Setting X = and observing^^'^^ that 

Jo r(f + £ + 



one obtains 

s ax^ r(£)r(a 
.^^ r(f + f + i) 



Finally, by inverting this formula one gets the Fermi energy Ep- □ 

THEOREM 2.2 Let us consider an ideal Bose gas in a power-law isotropic 
potential U{r) = A r" with r = |r| = (Z^^i a;^)^''^. The finite temperature 
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non-condensed momentum distribution is given by 



1 r(^ + i) / 1 \ " / 



The Bose transition temperature Tb reads 



B 



r(£ + i)c(f + f) 

and i/ie condensed fraction is 



D.D 



D.D 

No , /T\ - 



— = 1 - — 

where N is the number of Bosons in the gas. 

Proof. The finite temperature momentum distribution can be found by 
following the procedure used in the proof of the previous theorem: from 
the Eq. (6) and by integrating over space coordinates the Eq. (3) (but 
with the sign — ). It follows that one must use the Bose function go^z) 
with z — e^'''^~^\ At the BEC transition temperature Tb, the chemical 
potential // is zero and at // = the number N of particles can be analytically 
determined. One has 



- {2n^rv{^ + i)v{^ + i)\(3A) tii^^J^ 

where the Bose function has been written as a power series (see Remark 1). 
Setting X = iP-^ and observing^^'^^ that 



L 



dx x^e-^ = r — + 1 



2 
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and also 



=0 « 

one obtains 



^ 1 _.fD D\ 



2hV \Aj r(f + i) 



By inverting the function N — N{T) one finds the transition temperature Tb- 
Below Tb, a macroscopic number A^o of particle occupies the single-particle 
ground-state of the system. It follows that the previous equation gives the 
number N — Nq of non-condensed particles and the condensed fraction is 

No/N = 1 - (r/rB)^/2+^/". □ 

This last theorem generalizes the BEC results obtained with D = 3 by 
Bagnato, Pritchard and Kleppner.^^ 

It is important to observe that from the two previous theorems one easily 
derives the thermodynamic properties of quantum gases in harmonic traps 
and in a rigid box. In fact, by setting n = 2 and A = mu'^r'^/2 one gets the 
formulas for the Bose and Fermi gases in a harmonic trap (in the case of a 
anisotropic harmonic potential, cu is the geometric average of the frequencies 
of the trap). The results for a rigid box are instead obtained by letting 
^ —>■ 0, where the density of particles per unit length is given by N/Q^ and 

= Dtt^ /T(^ + 1) is the volume of the D-dimensional unit sphere. 

Finally, one notes that in the formula of the BEC transition temperature 
Tb it appears the function C(y -|- ^). Because ({x) < oo for x > 1 but 
= 00,-'^''''^^ one easily deduces the following corollary. 
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COROLLARY 1. Let us consider an ideal Bose gas in a power-law 
isotropic potential U{r) — A r"" with r = |r| = {J2iLi x^Y^'^ . BEC is pos- 
sible if and only if the following condition is satisfied 

D D ^ 
2 n 

where D is the space dimension and n is the exponent of the confining power- 
law potential. □ 

This is a remarkable inequality. For example, lox D — 2 one finds the 
familiar result that there is no BEC in a homogenous gas — 0) but BEC 
is possible in a harmonic trap (n = 2). Moreover, one obtains that for D = 1 
BEC is possible with 1 < n < 2. 



4 Conclusions 

By using the grand canonical ensemble of statistical mechanics and the semi- 
classical approximation, we have derived some thermodynamic properties 
of ideal quantum gases in a generic isotropic power-law confining external 
potential. We have calculated the density of states, spatial and momen- 
tum distributions and obtained analj^ical formulas for the Fermi energy, the 
BEC transition temperature and the Bose condensed fraction. Note that 
nowadays the spatial and momentum density profiles are quantities easily 
experimentally measured. We have also shown that BEC in an isotropic 
power-law potential is possible if and only \l D/2 -\- D/n, where D is the 
space dimension and n is the exponent of the confining power-law potential. 
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The present investigation is the starting point for future analyses of inter- 
acting quantum gases in D-dimensional space and generic trapping potentiaL 
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